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We study electronic structures at an interface between a topological insulator and a ferromag- 
netic insulator by using three-dimensional two-band model. In usual ferromagnetic insulators, the 
exchange potential is much larger than the bulk gap size in the topological insulators and electronic 
structures are asymmetric with respect to the fermi level. In such situation, we show that unusual 
metallic states appear under the magnetic moment pointing the perpendicular direction to the junc- 
tion plane, which cannot be described by the two-dimensional effective model around the Dirac 
point. When the magnetic moment is in the parallel direction to the plane, the number of Dirac 
cones becomes even integers. The conclusions obtained in analytical calculations are confirmed by 
numerical simulations on tight-binding lattice. 
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I. INTRODUCTION 

Physics of a metallic state on a surface of a three- 
dimensional (3D) topological insulator (TI)i~— is un- 
doubtedly a hot issue these days£~— . Intrinsic phenom- 
ena originated from the topological nature of the insu- 
lating state would open a novel field of condensed mat- 
ter physics. In particular, the metallic surface state 
shows interesting features when the TI is attached to an- 
other materials with gapped excitation spectra such as 
superconductors^ — and ferromagnetic insulators^—. 
The existence of Majorana fermion has been discussed 
in hybrid structures of such materials A 1 1 1 1 7 ~ — 

The surface metallic state has a linear dispersion so 
called two-dimensional Dirac cone. The upper and lower 
corns touch at a point in Brillouin zone, so called Dirac 
point. The 3D TPs can be classified into the strong and 
week TI in terms of the number of Dirac points ^ 2Q i 21 . 
Namely the metallic state is protected from the impu- 
rity scattering for odd number Dirac cone a 2Q i 22 , whereas 
it disappears for even number Dirac cones. To discuss 
intrinsic phenomena of TPs, it is necessary to tune the 
Fermi level near the Dirac point, which is possible in ex- 
periments by chemical dopin g 23 ! 24 . 

At the interface of a TI and a ferromagnetic insulator 
(FI), the metallic state is drastically modified depend- 
ing on the direction of magnetic moment 11 . The metallic 
state becomes insulating in the presence of magnetic mo- 
ment perpendicular to the interface plane. On the other 
hand, it remains metallic in the presence of magnetic 
moment parallel to the interface. The parallel magnetic 
moments only shift the Dirac point from the Y point in 
the Brillouin zone to another points there. Such con- 
clusions have been obtained by analyzing effective theo- 
retical model around the Dirac point, where the surface 
state is described by the two-dimensional Dirac Hamilto- 
nian under the small exchange potential due to the mag- 
netic moment. However it is unclear if these conclusions 
are still valid or not in real TI/FI junctions because the 



exchange potential of FI is much larger than the gap of 
TI. 

In this paper, we study electronic states at the interface 
of FI/TI junction by using three-dimensional two-band 
model. We show that asymmetry of the band structure 
in FI with respect to the Fermi level separates the dis- 
persion of interface state from bulk band in whole Bril- 
louin zone. This suggests that the effective theory around 
the Dirac point is no longer valid in real TI/FI junc- 
tions. A metallic interface state appears even when the 
magnetic moment in FI is perpendicular to the junction 
plane. When the magnetic moment in FI is parallel to the 
interface plane, number of Dirac points should be even 
number in whole Brillouin zone. In addition to large 
asymmetry of band structure in FI, breaking down the 
time-reversal symmetry and a basic feature of Brillouin 
zone also play important roles in there electric properties 
of the interface state. The conclusions obtained by ana- 
lytical calculation are confirmed by numerical simulation 
on three-dimensional two-band tight-binding lattice. Ob- 
tained results would be important not only in the basic 
physics but also in the view of potential application. 

This paper is organized as follows. In Sec. II, we 
summarize electric property at a TI/FI junction inter- 
face based on the effective Hamiltonian around the Dirac 
point. At the same time, we discuss the limits of the 
effective theory. In Sec. Ill, we analytically study effects 
of large band asymmetry and large magnetic moment of 
FI on the interface electric states. In Sec. IV, the conclu- 
sions based on the analytical results are checked by the 
numerical simulation on three-dimensional tight-binding 
model. The conclusion is given in Sec. V. 



II. EFFECTIVE THEORY AROUND THE 
DIRAC POINT 

We firstly summarize the features of the interface state 
which have been discussed by using effective Hamiltonian 



2 



around the Dirac point in two- dimension^. The effective 
Hamiltonian in two-dimension is derived from the three- 
dimensional electric states of a TI described by 



H 



Aso d(k) ■ s 
d(k) • s —Aso 

A=M -^B a kJ, 



(1) 
(2) 



where Mo and B a for a = 1 — 3 are band parameters. 
The unit matrix in spin space is denoted by so and s a for 
a = 1 — 3 are the Pauli matrices. The spin-orbit coupling 
is symbolically expressed by d(k) which satisfies 



d(-fc) = -d(fc). 



(3) 



The surface state on the TI is approximately described 
by the effective Hamiltonian in two-dimension, 



h S nv(k X: k y ) = v F D • 8 - fJL, 



(4) 



where vp is the Fermi velocity. In what follows, we im- 
plicitly consider E^Se^. However the arguments be- 
low are valid for all TPs. For Bi2Se3, it is shown that 
D = (—k y ,k x )£±. The Dirac point is at (k Xl k y ) = (0,0) 
which we call Ao in this paper. The dispersion relation 
becomes Ek = vp\k\ — fi. The spin configurati on on 
the Fermi surface is schematically illustrated in Fig. l(a)| 
where we assume /i > and focus only on the upper 
Dirac cone. The direction of spin and that of momentum 
are locked to each other. Thus the spin direction flips 
abruptly at Ao when we tr ace th e electronic states along 
the line L as shown in Fig. 1(b) . Thus the Dirac point 



may be a kink for the spin polarization on a line passing 
through it. This fact limits the validity of the effective 
theory around the Dirac point. Namely it is impossible 
to extend the effective theory to electric states in whole 
Brillouin zone. Let us trace electronic states along the 
straight line between Ai = (jr, 0) and A[ = (— 7r, 0) in 
the upper Dirac cone. The states at Ai and that at A[ 
must be identical to each other because the two points 
are connected by a reciprocal vector. In other words, 
the topology of the Brillouin zone is the same as that 
of two-dimensional torus (T 2 =^ x S 1 ). Although the 
energy of the two states are equal to each other, the spin 
direction of the two states are opposite to each other. In 
the effective theory, Ai and A[ characterize the different 
electronic states. In real TPs, the effective theory usually 
works well because electric states on the Dirac cone is ab- 
sorbed into the bulk energy bands before \k\ reaching at 
the zone boundary. 

The interface state between a TI and a FI is also 
approximately described by the effective Hamiltonian 
around the Dirac point in two-dimension, 



h T iFi(k x ,k y ) = h sur (k x ,k y ) + M • s 



(5) 



where M is the exchange potential in FI. Effects of the 
FI on the interface state are considered only through M . 



It is easy to show that the magnetic moment perpendic- 
ular to the two-dimensional plane, M z , gives rise to a 
gap energy at the Ao. The magnetic moment parallel to 
the interface (M x ,M y ,0), on the other hand, shifts the 
Dirac point from Ao to (M^/^p, — M^/^p). In addition 
to this, the fermi level stays at the Dirac point even in 
the presence of (M x ,M y ,0)l^. The conclusions obtained 
by analyzing Eq. ([5]) seem to be valid for weak exchange 
potentials smaller than the gap size of TI. However, the 
typical gap size in TI is 100 meV— i 26 i 27 7 whereas the gap 
of FI is of the order of eV— . Thus the low energy 
electronic states around the gap of TI should be studied 
by using more realistic theoretical model. 



III. EFFECTS OF BAND ASYMMETRY AND 
LARGE MAGNETIC MOMENT OF TI 

Let us consider a TI in three-dimension under the ex- 
change potential due to the magnetic moment in a FI. 
The Hamiltonian reads 



H _ ( h s d(k) • 8 
\d(k) ■ s —hoso 

h =M - B x k 2 - B 2 (k 2 + k y 2 ), 
d(k) =(A 2 k x ,A 2 k y ,A 1 k z ), 



(6) 

(7) 
(8) 



where Mo, A±, A 2l B±, and B 2 are material parameters 2 ^. 
The wave number in x, y and z directions are denoted by 
k x , k y and k Zl respectively. The Hamiltonian in Eq. (j6| 
is decomposed into two parts 



H=H + H', 

„ _({M -B l k 2 z )so A x k z s 



H' 



Aik z s z -(M -B 1 k 2 z )s J ' 

-B 2 (k 2 + k y 2 )s A 2 (k x s x + k y Sy) 
A 2 (k X S X ~\- ky Sy ) B 2 (k x 2 + k y 2 )s 



(9) 
(10) 

• (11) 



To analyze interface electric state, we apply the trans- 
formation k z in in Ho, 

_ ({Mo + Bik 2 )so iA 1 ks z \ ( v 
" 1 iMMz -(Mo + B^so ' [iZ) 



In Fig |l(c) , schematic band structures of Europium 
chalcogenides are illustrated. The band structures are 
generally asymmetric with respect to the Fermi level in 
these materials, which we consider through two param- 
eters Mi and M 2 with Mi ^ M 2 as shown in Fig. l(c)| 
The horizontal line shows the Fermi energy of FI. The 
lowest band and the highest one are spin-splitting due 
to the exchange potential. We assume that the middle 
bands are spin-degenerate. We consider the large asym- 
metry of the band structures through the the exchange 
Hamiltonian 



Hm = 
Mi =M + ii, 



Ms a + /n m so 




M 2 = M-/i m , 



(13) 
(14) 
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FIG. 1: (a) The spin configuration of the Fermi surface, (b) The spin configuration on the line L. (c) The schematic band 
structure of a ferromagnetic insulator. The arrow in a band denotes spin direction and the horizontal line means the fermi 
energy. 



where a indicates the direction of the magnetic moment 
in FI and /i m represents the asymmetry in the band struc- 
ture. In these definition, M\ = M2 and Mi ^ M2 de- 
scribe the symmetric and asymmetric band structure, re- 
spectively. 



A. perpendicular magnetic moment to plane 

When the magnetic moment of FI is perpendicular to 
the junction plane, the exchange Hamiltonian for the sur- 
face sate is 



Urn = 



/Mi 











f 


-M 2 




















{0 












(15) 



In usual FI's, a relation Mi ^> Mq holds. The Hamilto- 
nian Ho+H m is decomposed into two 2x2 matrices whose 
eigenvalues are E { {Mi - M )/2 with M x = M 1 + M 
and M2 = — M2 + Mo. The eigenstates of can be ex- 
pressed by 



vi(ki) = 





(16) 



The coefficients and bi satisfy 
h = 



{Mi + M )/2 - Bik? 



(17) 



where D\ = A\ and D2 = —A\. This Hamiltonian is 
equivalent to that of the surface state of a TI facing to 
vacuum by substituting (Mi + Mo)/ 2 by Mo. 

The imaginary wavenumber tzf takes different forms 
depending on the sign of Mi + Mq. For Mi + Mq > 0, k 



has the similar form as it is in the TI/ vacuum surface, 

(18) 




2B 1 (M i + M ) 



AS 



The eigenstate in this case can be described by 
'A Mi \ 



C-e 



-K i 



(19) 



with C± being arbitrary constants. For Mj + Mq < 0, 
the wavenumber becomes 




2B 1 (M i + M ) 



± 1 



(20) 



The eigenstate is given by 



= C< 



1 \ e -«i' * 



-Di/Ai | e _ Ki 



(21) 



For Mi > 0, Mi + M > always holds. Thus n x takes 
Eq. JTSJ). On the other hand, Mi + M can be either 
positive or negative even for M2 > 0. 

We first analyze weak exchange potential satisfying 
M2 < 2Mo- The wave function of in this case is 

(S) = ( Di/ i Al ) (^exp[- K + 2 ]+C7-exp[- K - 2 ]). 

(22) 

with Cf being the normalization constant. For simplic- 
ity, in what follows, we drop z dependence from the wave 
function. There are only two independent wave function 
for M2 < 2Mo. The surface state is a superposition of 
?pi and ip2 which are defined by 




-1 


i 



(23) 
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The total Hamiltonian Ho + H' + can be repre- 
sented in these basis of ipi as, 



M 1 \ fH^ H[ 2 
-M 2 ) + ^ £T£ 2 

Mi iv F (k x — ik y ) 

-iv F (k x + ik y ) -M 2 

with = ^2- The energy of the surface state is 



_ M l -M 2 



(Mi + M 2 ) 2 



VF 2 k 2 



(24) 
(25) 

(26) 



with /c 



fey 



For weak exchange potential 



M2 < 2Mo, the exchange potential in the z direction 
causes the gap, which is consistent with the previous the- 
ories^. The asymmetry of the band structures gives a 
constant energy shift to the dispersion relation. 

Next we consider strong exchange potential satisfying 
M 2 > 2Mq. In this case, the straight forward calculation 
of the eigenfunction at the T point results in 



4'i 



1 

72 



( x \ 



i 

w 



j ^2 



1 

71 



f \ 
1 



W 



, ^3 



1 

71 



-1 


i 



(27) 



For convenience, we employ an another basis as follows, 




V4 



'03 



(28) 



The total Hamiltonian i^o + H' + i^" m in this representa- 
tion reads, 

Mi -iv F (k x - iky) v F (k x - ik y Y 

H = I ^ F (^ + iky) -M 2 - 5 2 /c 2 



v F (k x + zfcy) 







-M 2 + B 2 k 2 



(29) 



with «f = A2/V2. The energy dispersion can be derived 
from the eigen equation, 

x 3 - 2Mx 2 - (B 2 2 k 4 + 2v F 2 k 2 )x + 2MB 2 2 k i = 0, 

(30) 

with x = E + M 2 . At the vicinity of T-point, x(k) is 
approximately given by 



x(k) = ao + ai/c 2 + a2& 4 . 



Here ao can be obtained easily by putting k = 0. We 
obtain two values as 



a 0, 2M. 



(32) 



For ao =0, a^ can be derived by putting the coefficients 
of k 4 and k 6 terms in Eq. [30] to be zero. Since M ^> Mo, 
a\ and a2 have simple expression 



ai 



v F 2 / v F * 

^M ±52 V 1 + i^W 



2M 



±£ 2 



(33) 

(34) 
(35) 



Then the energy dispersions are approximately given by 
E(k) = -M 2 ± B 2 fc 2 T J^2 B ^ 4 - ( 36 ) 



In the same way, we also obtain 



E(k) = Mi 



M 2M 3 ' 



(37) 



for ao = 2M. In both ao = and 2M, the coefficient 
of k 2 and that of k 4 have opposite sign to each other. 
In addition, we can also predict that two minima of the 
dispersion go across the fermi level and the interface be- 
comes metallic for M > 2Mq. 



B. parallel magnetic moment to plane 

When the magnetic moment of FI is parallel to the 
junction plane, the Hamiltonian of the surface sate at 
T-point is Hq + H m with 



Ms x + fi m S0 





(38) 



Here we assume that the magnetic moment is in the x 
direction. This does not loose the generality of argument 
below because the Hamiltonian is rotationally invariant 
in momentum space. Applying an unitary transforma- 
tion, we obtain 



Ui(H + H m )U = 
'(M + £ift 2 )s 

2<-^4i KjS x 

Mi 

M 2j 



M 



-(M + Bik 2 )so 



M 
with 



U 







_ ((so 



is v )/\/2 
(s 





is y )/V2^ 



(39) 
(40) 

(41) 



(31) The eigenvectors can be expressed by 



4h 




1p2 




(42) 
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FIG. 2: (a) The simple Brillouin zone and its TRI points A*. In this figure, A2 and A2 are same TRI points under translational 
operation of a reciprocal lattice vector, (b) The spin configuration on the line L with a single Dirac points is illustrated. 



The elements satisfy 



%A\K 



Mi + B 1 k 2 



(43) 



where M x = M +M and M 2 = M Q -M. The eigenvalues 
and eigenvectors can be calculated in the same way with 
the previous subsection. 

When the exchange potential is weak M 2 < 2Mq, the 
eigenvectors are The eigenvalues E\ are 



Ei = 



Mi _ _M 2 

2 ' 2 2 ' 

Corresponding vectors are given by 



1 



^ 1 = 71 



t~ x \ 






-1 

i 



(44) 



(45) 



The total Hamiltonian H = Hq + U^H'U becomes 



H 



Mi - M 2 



so 



M - v F k y 



-iv F k x 



(46) 



where vf = ^2 and 2M = Mi + M 2 . The energy disper- 
sion is given by 



E 



Mi - M 2 



±v F dk x 2 + (ky-M/v F )*. (47) 



The Dirac point moves from the F point to (0, M), which 
is consistent with the effective theory in around the Dirac 
point. The asymmetry of the band structures, however, 
shifts the fermi level from the Dirac point. 

When the exchange potential is sufficiently large sat- 
isfying M 2 > 2Mq, the basis of the surface state become 




^2 



^3 



(48) 



The total Hamiltonian H 
results in 



Hq + H' + H m in this basis 



H 



'Ml - V F ky 

iv F k x /V2 
-v F k x /V2 



-iv F k x /V2 
-M 2 + B 2 k 2 
—iv F k y 



-v F k x /V2 

iv F ky 
-M 2 - B 2 k 2 



(49) 



with vp = A 2 . 

By analyzing Eq. flj5| in detail, we can conclude that 

(a) there are two Dirac cones in the whole Brillouin zone, 

(b) the asymmetry of band structure in FI with respect 
to the fermi level may causes the separation of the inter- 
face state from the bulk band, and (c) three branches of 
surface states appear in the gap of TI. These conclusions 
can be confirmed in a simple case where we consider the 
dispersion relation along a line satisfying k x = 0. At 
k x = 0, three dispersion branches appear at the interface 

E 1 =M 1 -v F k y , (50) 

E 2 = - M 2 - v F ky^l + (B 2 2 /v F 2 )ky 2 , (51) 

E 3 = - M 2 + v F ky^l + (B 2 2 /v F 2 )ky 2 . (52) 

Near the T-point, two branches E\ and E 2 are almost 
parallel to each other. The remaining branch E3 goes 
across E\ and E 2 . Therefore there are two Dirac points. 
For k' y = k y - (Mi + M 2 )/2v F , Eq O and O can be 
represented by 



E = 



Mi - M 2 



±V F ky, 



(53) 



where higher order terms for k 3 in Eq (f52|) are ignored. 
The first term implies a asymmetry of the band structure 
of FI. 

As we have discussed above, the asymmetry of band 
structure in FI removes the dispersion of the interface 
state from the bulk band. This causes more drastic 
modification of interface state in the presence of mag- 
netic moment parallel to the interface plane. When 
M = (M^, 0,0), the magnetic moment shifts the Dirac 
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point f rom A q in the Brillouin zone to a point D as shown 
in Fig. 2(a) Let us consider the spin configuration along 
a line which satisfies D \\ M (Eq. [4j) and passes through 
the Dirac point D. For M = (M x ,0,0), the line corre- 
sponds to the stra ight line L connecting A2 and A 2 as 
shown in Fig. |2(a)| We note two key features of spin di- 
rection along the line: (i) A2 and A 2 are identical point 
to each other and (ii) the spin direction flips at D. If the 
number of the Dirac point is one, spin direction at A2 
and A 2 would be opposite to each other. This statement, 
however, contradict to (i). Therefore the number of Dirac 
point must be an even integer on A2 — A 2 . Since D is a 
Dirac point, at least one extra Dirac point is necessary 
on A 2 -A 2 (Fig|2(b)1). 

This conclusion above can be obtained in more gen- 
eral argument. The Dirac point can be regarded as the 
magnetic monopole in the momentum space. The Gauss 
integration in the first Brillouin zone becomes finite in the 
presence of the single monopole. This integration should 



coincide with the path integration of D(k) along the zone 
boundary. However the integration along the boundary 
vanishes be cause of the relation D(—k) = —D(k). Thus 
there must be extra monopoles in the Brillouin zone. Ac- 
cording to this argument, the number of the Dirac points 
must be even number in the Brillouin zone. In Eq. ([52]) . 
two Dirac points are expected in the present situation. 

The conclusions obtained by the analytical calculation 
are confirmed by numerical simulation in the next sec- 
tion. 



IV. NUMERICAL RESULTS IN 3D 



Let us consider a junction of TI and FI on three- 
dimensional tight-binding lattice as shown in Fig. |3(a)[ 
We describe the TI by using the two-band model as 



^=EE[ 2 L',i> g L',2. 



3d' k 



£tiso A ■ s 
A • s -^tiso 



Ck 



Ck,j,v,\. 

£ti =(Mo — 2b\ + 262 cos(k x a) + 262 cos(k y a) — 4b 2 — I^Ti)5j,j' 
A =(a 2 k x 5 j j,a 2 ky5 j j,-ia 1 (5 j j>+ 1 - 6jj>-i)), 



h(5 : 



jd'+l -r °jd 



(54) 



(55) 

(56) 
(57) 



where c J k ■ (ckj,u,s) is the creation (annihilation) operator of an electron with spin s, belonging to the band v = 1 — 2, 
having two-dimensional wave vector k = (k x , k y ), and at a lattice site j < in the z direction. We used the periodic 
boundary condition in the xy plane. 
In the same way, we describe the FI by 



od' fe 



(-£fi + E 2 )s + M • s 



c kdA 



(,fi =(2tcos(k x ) + 2tcos(k y ) — St — jJLFi)Sj,j' + + ^j/_i), 

for j > 0. At the interface (j = 0), TI and FI are connected by, 



Hl= Yl [4,0,1 '4,o 



(£/ + £1/2)50 A'-s 

A'-s (-C/ + E 2 /2)s Q + M-s/2 



Cfc,0,l 
Cfe,0,2 



2£,i =M - 26i + 2(b 2 + t) cos(k x ) + 2(b 2 + t) cos(k y ) - 4b 2 -St- fi TI - fi FI , 
2A! =(a 2 k x ,a 2 k y ,0). 

I 



(58) 
(59) 

(60) 

(61) 
(62) 



The hard wall boundary condition along with z-axis is 
employed. The parameters in this calculation take values 
of Bi 2 Se 3 : a x = 7.86M /a, a 2 = 14.6M /a, b x = 3.57 x 
10Mo/a 2 , and b 2 = 2.02 x 10 2 M /a 2 in TI side 26 . The 
lattice constant a is about 5 [A]. In FI, we assume 6pi = 
10" 2 6i and E 1 = -E 2 = -M/2. The total lattice size in 
the z direction is 200 sites, where TI and FI occupy 150 



and 50 sites, respect ively. A schematic band picture of 
a FI is shown in Fig |3(b)| Electronic structure becomes 
asymmetric with respect to the Fermi level. 

We first show the dispersion relations of the interface 
states rather large energy range for magnetic moment 
perpendicular to the interface (Fig. |4(a)[ ) and for mag- 
netic moment parallel to the interface (Fig. |4(b)[ ), where 
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FIG. 3: (a): TI/FI junction on the three-dimensional tight-binding lattice. The interface is at j = 0. (b): The schematic band 
structure of a FI. The arrow in a band denotes spin direction and the horizontal line labeled by cf is the fermi energy. 
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ky [7r/a] 




ky [n/a] 



FIG. 4: The global pictures of band structures are showed for a perpendicular (a) and a parallel (b) magnetic moment. 



the dispersion is calculated along k x = 0, fjbpi = —Mo/2, 
M 5M in (a) and M 2M in (b). The wave function 
of the interface states behaves like e J / JO for j < in TI 
with jo being the inverse of localizing length. In the Fig- 
ures, we also show the bulk band in TI. As we discussed 
in Sec. Ill, the upper dispersion in (a) is clearly sepa- 
rated from the bulk band of TI in whole Brillouin zone 
because of the band asymmetry in FI. The dispersions of 
the interface states for the magnetic moment parallel to 
the interface have rather complicated structure as shown 
in (b). We note that upper dispersion branche is well 
separated from the bulk band for \k y \ > 0.3. We zoom 
up the dispersion relations near the T point and discuss 
their features in the next figures. 

In Fig. 5(a)| we show the dispersion relation of the 
interface states along k x = for the magnetic mo- 
ment perpendicular to the interface. Here we assume 
fiFi = —Mo/2 and show the results for M = 2M and 
5Mq. When the magnetic moment is relatively small 
at M = 2Mo, the Dirac cone disappears as predicted 
by the effective theory around the Dirac point. When 



we increase the exchange potential at M = 5 Mo, how- 
ever, the dispersion of the interface state behaves like 
ek ~ olq — (12k 2 + o^/c 4 . As a result, the interface state 
become metallic. Features of the metallic also depends 
on the fermi level in the FI. The dispersion relation in 
Fig. |5(b)| show that the number of fermi surface is one 
for jipi = Mo j / 2, whereas for \ifi — — Mo/2 two fermi 
surface appears. These numerical results are consistent 
with analytical one's in Sec. III. 



Next we look into the interface states at TI/FI junction 
in the presence of the mag netic moment parallel to the 
junction plane. Figure 5(c)| shows the dispersion relation 
along k x = for M || x, where fiFi = and M = 2Mo. 
There are two Dirac cones in the Brillouin zone, which 
is consistent with the argument in Sec. III. In Fig. |5(d)| 
we show the results at M = 2Mo for fiFi = —Mo/2 and 
Mo/2. The characteristic features of the interface states 
are insensitive to parameters such as fi^i and M. 



8 



1.5 
1 

0.5 


-0.5 
-1 
-1.5 

1.5 
1 

0.5 



o 

K5 



-0.5 
-1 
-1.5 



(a) MFI = -M /2 




M = 5M - 
M = 2M - 



-0.5 



k y [ir/a] 



0.5 



(c) 




M = 2M ■ 



-0.5 0.5 

k y [it /a] 




fi F l = M /2 
= -Mq/2- 



0.5 



MFI = M /2 
/i FI = -M /2 



k y [n/a] 



0.5 



FIG. 5: The band structures of TI/FI junction with a perpendicular (a),(b) and parallel (c),(d) magnetic moment are plotted 
of the energy E versus the wave vector k y . The optical gap of FI is locked in 2Mo in (a) and (c). The magnitude of a magnetic 
moment of FI is 2Mo and 5 Mo. There are the surface band separated from the bulk band structure. The effect by shifting the 
Fermi energy in the optical gap is plotted in (b) and (d). 



V. CONCLUSION 

In this paper, we have studied electronic properties of 
interface state between a topological insulator (TI) and 
a ferromagnetic insulator (FI) by using two-band model 
in three-dimension in both analytically and numerically. 
The energy gap of FI is usually much larger than that 
of TI and the band structures in FI is asymmetric with 
respect to its fermi level. The dispersion branches of 
the interface state are separated from the bulk band in 
whole Brillouin zone due to the asymmetry of the band 
structures. When the magnetic moment is in the per- 
pendicular direction to the interface plane, the interface 
states become metallic. The number of fermi surfaces of 
such interface states depends on the material parameters. 
When the magnetic moment is in the parallel direction 



to the interface plane, metallic states always appear ir- 
respective of the amplitude of the exchange potential. 
The number of the Dirac point becomes even integers in 
whole Brillouin zone. Such drastic effects of the magnetic 
moment on interface states obtained in analytical calcu- 
lation have been confirmed by the numerical simulation 
on the tight-binding lattice. 
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